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Abstract 

We prove an approximate spectral theorem for non-self-adjoint opera- 
tors and investigate its applications to second order differential operators 
in the semi-classical limit. This leads to the construction of a twisted FBI 
transform. We also investigate the connections between pseudospectra and 
boundary conditions in the semi-classical limit. 
AMS subject classification numbers: 81Q20, 47Axx, 34Lxx. 



1 Introduction 

In the last ten years the theory of pseudospectra has developed rapidly, and has 
been shown to give substantial insights into the properties of non- self- adjoint (NSA) 
matrices and operators, [U El 13 UH 1221 12S|- In this paper we focus on its applica- 
tions to second order differential operators. This involves giving a new and more 
general definition of pseudospectra. Our first reason for extending the concept is 
that the standard definition does not provide any link with the geometry of phase 
space, which is of great importance in the theory of differential and pseudo differen- 
tial operators. By incorporating the connection into the definitions, we increase the 
conceptual clarity and facilitate the analysis of pseudospectra in those situations 
in which the semi-classical approximation is relevant. 

The second reason for concentrating on pseudo-eigenfunctions rather than pseu- 
dospectra is that the former are used in [Hj to provide a new method of solving 
evolution equations approximately. In several dimensions one could not hope to 
obtain sufficient pseudo-eigenfunctions by choosing just one for each point of the 
complex plane. Questions of spectral multiplicity arise just as they do for ordi- 
nary spectral theory, and indicate that a better parametrization is by points in the 
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classical phase space, not by complex numbers. We plan to use the results of this 
paper to extend those of to more general operators. 

The paper has three parts. In the first we prove an abstract approximate spectral 
theorem for NSA operators. We find a connection between this and quantization. 
The second part relates these ideas to the semi-classical analysis of differential 
operators via the semi-classical principal symbol of the operator and what we call 
interior pseudo-eigenvectors. Finally we introduce the concept of boundary pseudo- 
eigenvectors and describe how to construct them. We mention that ^24, contains 
results relating the boundary and interior pseudospectra of twisted Toeplitz oper- 
ators which are parallel to the ones which we obtain for differential operators. See 
|9j for related work on the wave equation. 

2 An Approximate Spectral Theorem 

In in] we have shown how to 'diagonalize' highly non- normal operators by using 
pseudospectra. The diagonalization is only approximate, but, in spite of this, 
it may be used to solve evolution equations efficiently for some quite singular 
infinitesimal generators. 

In this paper we formulate the underlying theorem at a general level, in order to 
make it accessible to a wider audience. All of the assumptions here are satisfied 
in the numerical examples discussed in jH], as we indicate in the next section. The 
ingredients are simple. We suppose that A is a bounded or closed, unbounded 
linear operator acting in a separable Hilbert space Ti. We also suppose that A 
is a multiplication operator acting in the space duo) where 1 < p < oo; for 

numerical calculations the simplest choice is p = 2, but p = 1 is more natural for 
some other purposes. We assume exphcitly that 



for all ip in the maximal subdomain of U'{VL)^ where the 'symbol' a : — C of the 
operator A is a measurable function and duo is a a-finite measure on VL. It is known 
that the spectrum of the operator A equals the essential range of a. We also assume 
that E : U'{VL) Tiis a bounded linear operator such that i?(Dom(A)) C Dom(A) 
and that 



(A^)(^) 



(T{uo)tlj{uo) 



AE-EA\\ < e 



(1) 



for a (preassigned, small) e > 0, in the sense that 





for all 4> e Dom(A). 
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Theorem 1 Let A be the generator of a one-parameter semigroup Tt acting on Ti 
and satisfying 

\\Tt\\ < Me^* (3) 

for all t > 0. Suppose also that 

Re{a{uj)) < 7 

for all u & fl. Then (OP implies 

\\TtE - Ee^^l < etMe^^ (4) 

for all t > 0. 



Proof Since the operators in Q are aU bounded it is sufficient to prove the estimate 
for aU (p G Dom(A). We then have 

r* d 



|Tt^0-Ee^Vll = II / -^m„,Ee^^</.) ds 

JO as 

< f \\Tt_s{AE - EA)e'''. 
Jo 

< [' \\Tt-s\\e\\e^'(l>\\ ds 
Jo 

< f Me^(*-^)ee^^||0|| ds 
Jo 

= e\\(f)\\tMe"". 



ds 



If y4 is a bounded normal operator then the spectral theorem states that one can 
find such a representation in which E is unitary, 6 = and the essential range of 
A equals the spectrum of A. The point of Theorem Q is that it may be applied to 
operators which are far from unitary and in situations in which the essential range 
of A is very different from the spectrum of A. The explanation of this relates to 
pseudospectral theory. 

One might try to develop an 'approximate functional calculus' based upon the 
above theorem. For example if Tt = e"^* is a contraction semigroup then under 
suitable conditions one can prove an analogue of Theorem for Ta^t = e"^""^^"* 
when < a < 1; see jHj. 

In order to compare Theorem Q with the results in [6| one needs to approximate E 
by an operator E' whose range is not contained in Dom(A). 

Corollary 2 If in addition to the previous assumptions one has \\E — E'\\ < e then 
\\TtE'(f)- E'e^'(f)\\ < e\\(f)\\{l + M + tM)e"" (5) 
for all G LP{Q) and all t > 0. 
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Proof This follows directly from 

+ \\T,{E-E')ct>\\ + \\{E-E')e^'ct>\\. 

The following modification of Theorem ^ assumes that one is given f & Ti and 
wishes to approximate Ttf . 

Corollary 3 If f ETi then under the conditions of Theorem^ 
\\Ttf - Ee'^Vll < 11/ - ^0||Me^* + £||0||tMe^* 
for all G L'P{9) and t > 0. 

Proof We have 

\\TJ - Ee^ll < ||Ti(/ - E(P)\\ + mE^ - i?e^V|| 
each of which is straightforward to estimate. 

The above results can only be useful if M, t and 7 are of order 1. There also has 
to exist (p such that ||/ — E(f)\\ and are both small. One cannot simply put 
(j) = E~^f, since E need not be surjective or invert ible. 

If p = 2, the standard way of solving this problem is to minimize the functional 

s{<f>) = \\f-E<j>r+6ur (6) 

for a suitable value of the regularization parameter 6 > 0; see ^2)- This is achieved 
in the numerical context by putting 

= E\(/©O) 

where E : L^{n) ^H® L'^i^) is defined by 

E<P = E<P® 51/2^. (7) 

We include the proof of the following well-known proposition for completeness. 

Proposition 4 If p = 2, the minimum of (0j is achieved for = Egf , where 

Fs = {E*E + 6I)-^E* (8) 

satisfies \\Fs\\ < 6^^^"^ . Moreover \\EFs\\ < 1 for all 5 > 0. One has 

hmEFsf = f (9) 

for all f eH if and only i/Ran(£') is dense in Ti. 
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Proof The first statement depends upon a routine variational calculation. For the 
second we observe that 

11(^*^ + 5/)"^^" < ab 

where 

a= \\{E*E + SI)-'/^\\ < r^/' 

and 

b^ = \\{E*E + 5I)-^/^E*\\'' 

= II {E*E + 61)-^'^ E\ E{E*E + 

< \\{E*E + 6I)-^^\E*E + 6I){E*E + 61)-^^^ 

= 1. 

This calculation also implies that 

IIEF^II = \\E{E*E + 6iy^^\{E*E + 6I)-^/^E*\\ 
= \\{E*E + 6I)-^^^E*f 
< 1. 

Since Kan{EFs) C Ran(i?), Q implies that Ran(£') is dense. If Ran(£') is dense 
then the uniform boundedness just proved implies that holds for all / G 7i 
if it holds whenever / = E(j) for some G In this case let P denote 

the orthogonal projection onto the closure of the range of E*E. Since Ker(£^) = 
Ker (£'*£'), we may assume without loss of generality that P(f) = 0. We have 

Mm EEs f = \imE(E*E + 6I)-^E*E(j) 
= EP(f) = E(l) = f 

by applying the spectral theorem to the non-negative self-adjoint operator E*E. 

Using Proposition|3]one may ensure that e\\(f)\\ is small by choosing 6 appropriately. 
Even if E has dense range, one cannot ensure that ||/ — E(f)\\ is small for some 
particular 6 > without further conditions. One has either to make the a priori 
assumption that / lies in some subspace of well-approximable vectors, or observe 
a posteriori for particular choices of / and S that the minimizing does indeed 
make this quantity small enough for the application intended. 

3 The Connection with Pseudospectra 

Given e > 0, the e-pseudospectrum of the closed operator A is defined by 
Spec,(A) = {z : \\Af - zf\\ < e\\f\\ for some / G Dom(A)}. 
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Pseudospectral ideas lie at the core of this paper, and we refer to 001311111221122] 

for background material on this subject. The following theorem is valid for all 
p G [1,cxd), but its main apphcation is for p = 1. Indeed we conjecture that if 
p = 2 the first condition on E can only hold if E is isometric. In the following 
theorem Pu denotes the operator of multiplication by the characteristic function 
of the set U, always assumed to be measurable. 

Theorem 5 Suppose that 1 < p < oo, \\EPij\\ = 1 for all subsets U of Q with 
positive measure, and \\AE — EA\\ < e. Then 

Spec(A) C Spec,(A). 

Proof Let {3 G Spec (A). We choose 5 > such that 

e' := \\AE - EK\\+5 <e 

and put 

U = {uje^: \X{lo)-P\ < 5}. 
If (j) has support in U then 

\\AEcj>-i3EcP\\ < \\iAE-EA)<p\\ + \\EiA<p-p<P)\\ 

< mi 

Therefore 

M{\\Af-f3f\\/\\f\\:0^fen} < mf{\\AE<j>-(3Em\\E<f>\\:0^<peLP{U)} 

< e'mf{\\4>\\/\\E^:0j^^eL^m 
= e' < e. 

This implies that /? G Spec£(y4). 

Theorem 6 If p = 1 then the conditions of Theorem [3 are equivalent to the fol- 
lowing statements. There exists > and a set N of zero measure, such that for 
each UJ G VL\N there is a unit vector e^^ G Dom(A) which depends measurably on 
oj and satisfies 

\\Ae^ - a^e^W < e' < e 
for all UJ G Q\N, where <j{uj) G C. 



Proof The passage from the assumptions of Theorem El to the statements of this 
theorem is justified by using ^21 Theorem VI. 8. 6]. 

If lies in the maximal domain of A then under the assumptions of Theorem IHl 

\\AE(P - EA<p\\ < f \<f){u)\ \\Ae^ - a^ej du < e'Ul (10) 
Jn 
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Hence \\AE — EA\\ < e' < e. The calculations involved would be easy to justify if 
one only had to deal with finite sums, or if A and A were bounded, but in general 
they use limiting processes to define the integrals. Commuting A and A with these 
limiting processes is justified by the following lemma. 

Lemma 7 Let A be a closed linear operator with domain in a Banach space B and 
range in a Hilbert space Ti. Let c> , fn & Dom{A), — /|| — 0, \\gn — g\\ — ^ 
and \\Afn ~ gn\\ < c for all n, then f G Dom(y4) and \\Af — g\\ < c. 

If Q has finite measure then L^(fi) is continuously embedded in L^{Q), and all 
of the theorems of Section 1 hold under the present hypotheses. In the numerical 
applications of |S] the space fl is taken to be the finite set {!,.., A^} and du is the 
counting measure. Given unit pseudo-eigenvectors e„ G 7i of A for 1 < < A^, we 
have 

N 

E0=^0„e„. (11) 

n=l 

There is no requirement that the vectors should be linearly independent, and indeed 
in some of the examples studied in they are taken from an overcomplete infinite 
sequence {e„}5^^. Equivalently the operator E need not be invertible, or may have 
a large condition number. 



4 Quantization 

In this section we make some general comments about the relationship between 
our previous results and the notion of quantization. 

Let f2 be a second countable locally compact Hausdorff space, and let du be a 
regular Borel measure on Q with support equal to Q. Let 7i be a separable Hilbert 
space and let e : — > 7i be a continuous function. We define E : Cc{^) H hj 

E(f) = / (j){uj)ec^ duj. 

The following are well-known and elementary. 

Lemma 8 The operator E extends to a bounded linear operator Ei : L^{Q, duj) 
Ti if and only if u — >■ ||e(^|| is a bounded function, in which case 

\\Ei\\ = sup{||e^|| : u G fi}. 

The operator E extends to a bounded linear operator E2 : du) ^ H if and 

only if 

[ \{f,e^)\'duj<cW (12) 
Jo. 
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for some c > and all / G 7Y, m which case II-E2II is the smallest such constant c. 
The operator E* : Ti C{Q) is an isometry from Ti into L'^{Q, du) if and only if 



for all fen. 

Families of vectors {e^}^^fi satisfying (fT^ are also called continuous resultions of 
the identity and have played an important part in group representation theory and 
quantum mechanics for many decades. For their connection with coherent state 
theory and the Bargman transform see |21 ch. 8] and fTfi ch. 3]. If fll3|l holds then 
E\f = E* f for all / G 7Y. but this is not the case under the assumption (|T^ . 
which is more relevant to this paper. 

Given a function / G Cc(fi) we define the multiplication operator Mf by Mfcf) = fip 
where G U'{Vt) for some p. We define the quantization of the function / to be 
the operator Q{f) = EMjE* on Ti. We may also write 



where Paip = {tp,a)a] see, for example, j21 sect. 8.5]. The following lemma is also 
standard. 

Lemma 9 // / > then Q{f) > 0. If E is bounded from L^{n) to H then Q 
extends to a bounded linear operator from L^{Q) to the space T{li.) of trace class 
operators on Ti. If E is bounded from L'^{Q) to Ti then Q extends to a bounded 
linear operator from L°°{VL) to the space C{T-C) of bounded operators on Ti. Given 
/ f73|) . or equivalently EE* = 1, we have Q(l) = 1. 

In quantum theory it is commonplace to refer not to the operator Q but to the 
positive-operator-valued measure A{U) := EM^^E* where xu is the characteristic 
function of the measurable set U of Q. The formula 



implements a one-one correspondence between the two definitions; see [21 Lemma 
3.1.2]. If EE* = 1 then A{Q) = 1 and A{-) is called a generalized observable; for 
a systematic study of POV measures and their relation to coherent states see j21 
Ch. 3] or . See for more recent references and a connection with subnormal 
operators. 

The difference between this method of quantization and the approach of this paper 
is now clear. Instead of studying Q{f) = EMfE*, we would like to study S{f) = 




(13) 




Qif)= [ f{cu)A{du;) 
Jn 
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EMfE^^. If this were possible / — » S{f) would be an algebra homomorphism 
from L°°{Q) to CiTi). Since E is not invertible in general we compromise by 
studying EMfFs, where the regularized inverse Fs is given by (jH} and 5 > is 
chosen small enough to yield numerically valuable results but not so small that the 
computational algorithms become unreliable. 

The operator E which we have considered above has much in common with the 
Fourier-Bros-Iagolnitzer (FBI) transform as defined in jTTl ch. 3]. See also Ch.3], 
where the connection with the Wigner distribution and applications to quantum 
theory are explained. In Section |H1 we define a distorted FBI transform; the distor- 
tions are introduced to adapt the transform to a given differential operator, and 
involve replacing the Gaussian states used in the definition of the FBI transform 
by pseudo-eigenf unctions of the operator. 

5 The Connection with Semi-classical Analysis 

Before describing the connection of the above ideas with semi-classical analysis, 
we generalize the notion of pseudospectra. Following [71 HH Uni 1211 , define the 
(generalized) pseudospectra of a family of closed operators {A^}^^^q acting from 
dense domains Dom(y4^) in a Banach space B to another Banach space C to be the 
sets 

Spec,(A) = {lu : \\A^f\\ < e\\f\\ for some / E Dom(A^)} 
where e > 0. We have 

Spec,(A) U Spec,(v4*) = Spec(A) U {u : pj^ > e'^} 

where Spec (A) is defined to be the set of uj for which is not invertible. If 
dim(i3) = dim(C) < oo then 

Spec(A) C Spec,(A) = Spec,(A*) 

for all e > 0. If dim(i3) < dim(C) < oo then 

Spec,(A*) = Spec(A) = Q 

for all e > 0. The proof of the following lemma may be found in [T^ . 

Lemma 10 One has u G SpeCg(A) if and only if there exists a bounded operator 
D : B ^ C such that \\D\\ < e and 

KeT{A{uj) + D) ^ {0}. 
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Given a differential or pseudodifferential operator Lh with domain C^(X), where 
X is a region in and /i > 0, we define the operator family 

by 

Ah,u,if = {Q'f - u^f, P,f - 0/, Lhf - ct{u, Of) (14) 

where {Q^ f){x) = f{x) and {Pjf){x) = —ihdjf{x). In these equations we assume 
that M e X, ^ e R^, 1 < j < iV and (t(m, ^) is the semiclassical principal symbol 
of the operator L^, as defined below. It follows directly from the definitions that 
\\Ah,u,if\\ < implies 

< ^11/11 
\\Pjf-Of\\ < ^11/11 
\\L,f-a{x,Of\\ < e\\f\\ 

where 1 < j < X. It is known that the pseudospectra converge to fill a certain set 
(t(A) if /i — > and e ^ simultaneously at suitable rates; see Section 7 for details. 
Even in one space dimension a point in cr(A) may be the image of more than one 
point in A, so cr(A) may have hidden structure as a subset of C. This observation 
applies with less precision to the numerically determined pseudospectra for fixed 
h> and e > 0. 

The extension of the above ideas to a manifold X needs some care, since the full 
symbol ah{u,C,) is not an invariant object in general. It is shown in (20] that one 
can resolve these problems if the manifold is provided with a linear connection, as 
happens if it is Riemannian. The symbol ah{u,C,) is then definable as a function on 
the cotangent bundle T*X and A is a certain subset of T*X. We do not actually 
need the full symbol for our problem: its semiclassical limit is sufficient. The 
semiclassical principal symbol is given by 

n— >0 

and is an invariant quantity, function on the cotangent bundle T*X it 

does not depend on the choice of local coordinates. 

The following alternative definition of the semiclassical principal symbol of Lh 
makes its invariant character clear. Suppose that u E X and ^ is a cotangent 
vector at u. Let / be any smooth function on X such that df{u) = ^. Then 
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6 The Semiclassical Spectrum 



The theory which we shall describe can be developed at several levels of generality, 
and in this section we consider only second order differential operators acting on 

Given /i > 0, let denote the operator 

(Um^) = -h'ai\x)dj,kf{x) - ihl/^{x)djf{x) + Ch{x)f{x) 

acting on functions / : — > C, where a, b, c are sufficiently regular functions 
whose values are respectively matrices, vectors and scalars with complex-valued 
entries, and we use the standard summation convention. Under conditions which 
we shall impose the domain of will contain C^(R^). All considerations in this 
paper are local, so no growth bounds at infinity on the coefficients are needed. We 
allow the coefficients to be /i-dependent so that the class of differential operators 
is invariant under local changes of coordinates. The semiclassical principal symbol 
of this operator is the complex- valued function 

a{u, = a^o''{u)ijik + &o(«)Ci + Co(m) (15) 
in which we take ^ to be real vectors in R-^. 

Given {u,^) e R-^ x R-^ we are interested in finding locahzed approximate eigen- 
functions for the operator L^. We require that they become asymptotically exact 
as /i ^ 0. 

Our first theorem provides the motivation for defining the semi-classical spectrum 
of Lh to be the set a(R^ x R^). 

Theorem 11 Suppose that o^j^{x), and Ch{x) are all locally Lipschitz con- 

tinuous in both X G R^ and h G [0,1]. Then for every u G R''^, ^ G R''^ and 
h G (0, 1] there exists fh G ^^(R^) such that 



IIMI2 = 00 (16) 

WQ'fh-u^fhh = 0{h'/') (17) 

\\Pjfh-Cjfhh = 0(/iV2) (18) 

Lhh-<j{x,0fh\\2 = 0{h'/') (19) 



ash^O, for all l<j<N. 

Proof Let be a non-negative C°° function on R^ which equals 1 if \x\ < 1 and 
if \x\ > 2. Given {u, ^ e R^ x R^, /i > and a = 1/2 define 

fh{x) = /l-^°/2gift-ie.x^(^-a(^ _ 
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The first three statements of the theorem are routine verifications performed by 
the same method as follows. 

We verify (|T!Hl by using the expansion 

Lhfh - OA = g + ri + r2 + rs + r4 

where 

g(x) = {a^h'^i^) - O'o''{u)}^j^kfh{x) 

+{bi{x) - KiuM^Mx) + {ch{x) - co{u)}Mx) 

and 

n = -ih^~''-^''/'^a^'\x)^je'''''^-''(f)k{h~"{x - u)) 
r2 = -th^-''-''''/^a''\x)^ke''''"^-''(j)j{h"'{x-u)) 

r4 = -z/ii-"-^"/26^(x)eje"'"'«->,-(/i-"(a;-u)). 

In these identities the subscripts on (p denote partial derivatives. The Lipschitz 
assumptions on the coefficients of L^, and the fact that the support of fh has 
diameter of order imply that 

\\9h = 0m 

as h ^ 0. We also have \\rj\\2 = 0{h^~°') for j = 1,2,4 and ||r3||2 = 0(/i^~^°). 
The overall error is minimized by putting a = 1/2. 

7 Constructing the Interior Pseudospectra 

The material in this section is based upon the fact that if the coefficients are 
sufficiently smooth then the estimate ()19|) can be greatly improved by a suitable 
choice of fh- In the language of Section 5 we replace (fT^ by 

Ak,u,J = {Q'f - u^f, Pjf - 0/, h-^LJ - Mu, Of}) (20) 

where n > 0. The size of n depends upon the smoothness of the coefficients, which 
for simplicity we assume to be C°°. The pseudospectral estimate < £:||/|| 

then implies 

WQ^f-u^fW < e 11/11 

\\P,f-Of\\ < e\\f\\ 
\\Lhf-a{x,Of\\ < h^eWfW 
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where 1 < j < A^. We repeat the calculations of [3 E] for a more general second 
order ordinary differential operator for completeness. The extension to pseudo- 
differential operators in higher dimensions, [HlEnj, cannot be formulated in exactly 
the same manner: there can be infinitely many different pseudo-eigenfunctions 
associated with a point in phase space, and the correct parametrization of these is 
not obvious. We assume that 

iLJ)ix) = -h^a{x)f\x) - thb{x)f'{x) + c{x)fix) 

so that the semiclassical principal symbol is 

a{u,0 = aiu)e + biu)^ + ciu). 

We assume ellipticity, in other words that a{x) ^ for all x G R. Given m, ^ G R, 
we put 

fiu + s) = h-'/\{s)exp{i;is)) (21) 

for all s G R, where x ^ satisfies x{s) = 1 if |s| < 6/2 and x{^) = if |s| > 6, 
and 6 > must be small enough; see the proof of Lemma [T21 We assume that 

n 

^P(s) = Yl h^i^mis) (22) 

m=— 1 

for some integer n> —1. This is a non-standard form of the JWKB expansion, and 
has the feature that the function / does not vanish within the interval of interest. 
A direct computation shows that 



(2n+2 \ 
E ^"'^^ / + Rem 
m=0 / 



(23) 



where Rem = 0{h°°) a.s h —>■ under the conditions which we impose below. Also 

0o(s) = -a{u + s){^L^{s)f -tb{u + s)tp'_^{s) + c{u + s) 
—a(u)^^ — b(u)^ — c{u). 

Assuming ellipticity, that is a{x) ^ for all x G R, the eikonal identity 0o = 
implies 

Jv=o [ 2a[u + V) ^ J 

where 

a(M)^2 b{u + vy c{u)-c{u + v) 

Wiu, 4, V) = — H -. + -— TTT H 7 r . 

a[u + V) a[u + v) Aa[u + vy a[u + v) 

We take the branch of the square root which equals ^ + b{u)/2a{u) at f = 0. The 
condition implies that da/d^ ^ and hence that w(m,^,0) is non-zero; this 
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implies that w{u,C,,v) ^ for all small enough v; and hence that the square root 
is uniquely determined for all such v by the requirement of continuity. 

Writing ip-i{s) = i^s + ks'^/2 + 0{s^) for some G C, we then obtain 

-ik{2a{u)^ + b{u)} + a\u)^'^ + h'{u)^ + c\u) = 0. 
The requirement that Re (k) < may be rewritten in the form 

and then in the form (m, ^) G i7 where 

n = {{u,O:Wu(T,}>0} (24) 

and 

r 1 ._ dcTi d(T2 

du du 

and (Ji = Re (cr), (T2 = Im (cr). In examples one may find that Vt is not connected. If 
it has components Aj then ij{VLj) may overlap. The multiplicity of a point z G a{VL) 
may be defined by 

If the coefficients of are smooth then for any choice of n one may choose ipoy ■■■y'4^n 
so that 01 = ... = 4>n+i = 0. This is achieved as follows. If 1 < m < n then 



->m+l 



It follows from (j24j) that 2ail)'_i + ife 7^ if s = 0, and hence that it is non-zero for 
all small enough s. If we define ipm by 

^™(«) = / — ^7 

JO 2aip_i + 

Then l^/^mls)! < Cm|s| and = for all small enough s. On making these 

choices we obtain a pseudo-eigenfunction /, depending on h,n,u and ^, for which 
Lhf - a{u, 0/ = 0(/i"+2) as /i ^ 0. 

The proof of Theorem IT^ below is facilitated by introducing the scale of spaces S'^, 
consisting of all functions which can be written as finite sums of functions of the 
form g(s) = h°'~^/^s^ p{s) exp{?/'(s)} where ip is given by (j22p . p G has support 
in [-5, 5], a G R, /3 G {0, 1, 2, ...} and 2a + /5 > 7. Putting 8^ = n^ea^^ we see 
that if, in addition to the above assumptions, p vanishes in some neighbourhood 
of 0, then g G £°°. 
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Lemma 12 If 6 > is small enough and g G S"^ then there exists c such that 



for allO<h<l. 



Proof It is sufficient to consider the case in which g is one of the terms of the 
form assumed in the definition of £"'. One may rewrite \h^^^~°'g(s)\'^ in the form 
s^'^G(s) exp{— where F{s) = — 2Re {ip_i{s)/ s"^) is a positive continuous 
function on [S, 5] if 5 > is small enough and G is a continuous function on [—6, 6]. 
By Laplace's method we have 

/' s^''Gis)exp{~h-'s^Fis)} ds ~ ch^^^+'^/^ 

duS h —> 0+, where 

G(0)r((2/3+l)/2) 
^ F(0)(2/3+i)/2 

The statement of the lemma follows immediately. 

Theorem 13 If the coefficients of Lh are C°° and (u,^) lies in the set Q defined 
by l{24\ ), then for every positive integer n there exist functions f G depending 
on h, n, u, ^ such that 

hm 11/11 = 00 (25) 

WQf-ufW = 0{h'/') (26) 

\\Pf-^f\\ = 0{h'/') (27) 

\\L,f-a{u,Of\\ = 0{h-^') (28) 

as h ^ 0. 



Proof We define / by (j2^ and observe that f E £^ . The asymptotic formula (j25|) 
follows by the method of proof of Lemma El We next observe that Qf ~ uf E 
so (j26j) follows from Lemma IT^ 

We have 

Pf - ^f = 1^1+ 1^2+ /is 

where 

fii = -^h-'/^{^PUis)-^Oxis)exp{iJ{s)}eS' 

n 
m=0 

/i3 = -th'/Ws)exp{ij{s)}eS°^. 
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Therefore Pf — ^/ G and (j2Z|) follows using Lemma IT^ 
Since (pm = for < m < n + 1 it follows from ()23|1 that 

/ 2n+2 \ 
\m=n+2 / 

This implies (f^Hj) by Lemma IT^ 

Note The orders of magnitude of the errors in ()26|) and ()27|) cannot both be reduced 
by a different choice of the function /, because of the uncertainty principle. 

The following lemma shows that one can approximate the pseudo-eigenfunction by 
a Gaussian expression. 

Lemma 14 We have 

\\f-9\\<ch'/' 

as h ^ 0, where 

giu + s) = h-^/^ e^^{h-^{iis + ks^/2)}. 
Proof Since g — X9 = 0{h°°) we have to estimate the norm of 

By virtue of the bound 

|g-a _ ^ _ ^|g-min(Re(a),Re(b)) 

this is dominated by the absolute value of 

^(s) = h-^/^{4j{s) - h-\i^s + ksy2)}x{s) exp{-h-^cs^} 

for some c > 0. In the following calculations we define in the same way as 
but with ip{s) replaced by —h~^cs^. We may write /i = /ii + /i2 where 

f^2{s) = h-'/^(j2h'^Ms)]x{s)exp{-h-'cs^} 

\m=0 J 

Since 

we have fj.i G S^. Since < Cm|s| for all s we also have /i2 G S^. The estimate 

of this lemma now follows by an obvious modification of Lemma ^1 
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8 A Semi-classical Transform 



We continue with the assumptions and notation of the last section. Theorem El 
provides the information needed for the apphcation of Theorem IHl We define the 
set Q in Theorem IHl by (j^^ and take a to be the semi-classical principal symbol 
(fT3|l of A. In numerical applications, one would, of course, have to restrict to a 
finite subset of fl, as described in jH]. 

We fix n and put eh,u,^ = fh,u,(/\\fh,u,^\\ where fh,u,^ = / is defined by The 
semiclassical integral transform E : L^{Q) -^^(R-) is then defined by 



{Eh(j))ix) = J^(j){u,^)eh,u,^{x) dud^ 



and has norm 1 by [iDi Theorem VI.8.6]. The functions fh,u,i{x) are very compli- 
cated for large n, and the following approximation may therefore be valuable. 

Theorem 15 Given h,u,^, let 

9h,u,i{x) = h-'/' exp{h-\iCix ~u) + K,^{x - uf/2)} (29) 

where 

da { da\ , ^ 

// (m,0 e n then gh,u,i e L\R). Define E'^ : L\n) L^{R) by 

{E'^c^){x) = f 0(n,Oe;,„,^(x) dud^ (31) 

where e'^,^^^ = gh,u,i/\\gh,u,^\\- Then \\E'f^\\ = 1 and 

hm\\Eh<P-E'f,<p\\=0 (32) 

for all (f) G L^{Q). If we replace Q by a compact subset U of Q then 

^m\\E,,-E',\\=0 (33) 



Proof We start by observing that Re (/c„^g) < if and only if (u,^) G Q, so 
gh,u,£, ^ -^^(R) under the same conditions. We have = 1 by |T01 Theorem 

VI.8.6]. 

Let {fln}'^!^=i be an increasing sequence of compact subsets of fl whose union equals 
Q. If we can prove that the restrictions E^ n and to L^(f2„) satisfy 

lim||E,,„-E;_J| = (34) 
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then and follow by standard procedures. 
In Lemma IT4l we proved that 

\\fKu,i-gh,u,d = 0{h^'^) 

for each (m, ^) G as /i — * 0. The dependence of the error upon m,^ and the 
coefficients of A was given explicitly, and implies that 

limsup{||/h,„,^ - 5f/,,„,5|| : ^ Vtn} = 0. 

Taking into account we deduce that 

limsup{||e;,,„,^ - e'f^^^J : G ^n} = 0. 

This implies (jHlI). 

Lemma 16 Let let E'^^ denote the restriction of E'f^ to the subset U of Q. If 
U, V are two compact subsets of Q which are spatially disjoint in the sense that 
(u,^) G U and {v,ri) G V implies u ^ v then the ranges of E'^jj and E'^^ y are 
uniformly asymptotically orthogonal in the sense that 

lim||(E;,,^)X,yll=0. 

The convergence is exponentially fast. 

Proof Let W be an open subset of R such that f/ C {W xK) and (Wx R) = 0. 
Let P be the orthogonal projection in i^^(R) whose range consists of all functions 
with support in W . Then 

WiKuTKyW < \\{Elu)\I-P)Ky\\ + \\{E',^vrPKy\\ 

< \\{I-P)El^\\ + \\PEly\\. 

We consider further only the first term on the RHS; the other is treated in a similar 
manner. If G L^{U) then 

Wil - P)E',<P\\ = \\J^{I~P)el^^^<Piu,Oduda 

< Jj{l-P)el^jmu,0\dud^ 

< snp{\\iI-P)el^J:iu,O^U} 
snp{\\iI-P)g,,ud--iu,O^U} 



< 



inf{||^;,,„,^|| : e t/} 
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The explicit expression ()29p for g and the compactness of U ensure that the final 
supremum converges to exponentially fast as /i — ^> while the final infimum 
converges to a positive limit. 

If we subdivide R into small intervals then the lemma implies that E'^ (or more 
exactly its restriction to any compact subregion of f2) acts asymptotically indepen- 
dently on subintervals which are not adjacent. If each interval is small enough we 
may approximate E'f^ in any subinterval by the operator with a frozen value of u. 

We conjecture that under suitable conditions on the coefficients of A the transforms 
Eh and E'^ are both bounded from L^(f2) to L^(R). As evidence for this we treat 
the case in which the variable u in fc^^^ is frozen at the value v. We also assume that 
A is a Schrodinger operator, so that its symbol is of the form (t(m,^) = + c{u). 
This implies that ky^^ = —1/nC, where k = 2/ic'{v). Assuming that k has positive 
real part, it is immediate that Re ky^^ < if and only if ^ > 0. We therefore put 
R2 ={(u,O:«eR,e>0}. 

We define the distorted FBI transform E^ : Cc(R^) -^^(R) by 

Eh(t> = h-'/^ [ 0(m, Oe/.,n,5 du de (35) 

where eh,u,^ = gh,u,^/\\9h,u,^\\ and 

gh,u,ii{x) = exp{ii^(x - u)/h - (x - uf/2hKQ. (36) 

Theorem 17 //Re (k) > and h > then the operator 0^5]) may be extended to 
a bounded operator from L^(R^) to L'^(R,) whose norm is bounded above uniformly 
as h ^ 0. 



Proof In this proof we write Cr to denote positive constants which depend only on 
K. We always take ^ to be positive. We have 

= J^exp{-Re{l/K){x-uY/h^} dx 

Therefore 

\\~gh,uA\ = C2h'/'e^\ 

We prove the boundedness of E^ rather than that of Eh- We have 

iE*hf)iu,0 = I K{u,^,h,x)f{x) dx (37) 
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where 

K{u,^,h,x) = h'^''^eh,u,(,{x) = (3h,il£,{u- x), 

7g(M) = exp{i^u/h — /2hK^}. 

We next take the Fourier transform JF of ()37|) in the u variable, noting that ^ is a 
unitary operator on L^(R^). This yields 

= \\k\\ 

where 
and 



7S 



We deduce that 

\\Elf\\<c,\\:Ff\\=c, 
for all / G /^^(R) if and only if 



Our task therefore, is to prove that the function 

/■oo 

= / h~'/^e'^ exp{-C6(e//^ - 5)^/^0 (38) 

JO 

is bounded on x R, provided Cg > 0. If s < then putting ^ = h^^^r] we obtain 

< F(/i,0) 

/■oo 

Jo 

= / r/^'^ exp{— C6?7'^} dr], 
Jo 

which is finite. If s > then putting = hsrj we obtain 

F{h,s) = G{hh^) (39) 

where ^ 

G(t) = I r/V2ii/2g-c6(„-i)^,*^^^ 
Jo 
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so we have to prove that G is bounded on (0, oo). We do this in stages. If < t < 1 
then 





because every term in the integrand is less than 1. If t > 1 then 
Jo 

roc 



lit>0 then 



4 







Cjt ^ < C7. 



1/2 Jl/2 

J —00 



-00 



Finally if t > then putting i] = (t we obtain 

;^i/2^i/2g-c6(,-i)2,t^^ < j^-'/H^I^-'^^'Ur^ 



00 



< / C^/^e-«^^dC 

JO 

= Cg. 

One cannot expect i?^ to be isometric, as is the case for the FBI transform, but we 
prove that this is asymptotically true in the semi-classical limit, up to a normalizing 
constant c, which could be evaluated explicitly. 

Theorem 18 There exists a positive constant c such that 

\iT^\\Elf\\=c\\f\\ 

ft— >L) 

for all f G L2(R). 

Proof In the proof of Theorem El we obtained the formula 

/oo 
F{h,s)\{Tf){s)\' ds 
-00 

where 

< F{h,s) < cii 
21 



for all /i > and s G R. By the dominated convergence theorem it suffices to 
prove that 

lira F(h,s) = C12 := / 77^'^ exp{— cer/^} di] 

for all s G R. We do this for s > 0, noting that the cases s = and s < are 
similar. By it suffices to prove that lim(_^o+ G{t) = c^. As t ^ 0+ we have 

G{t) = r r^^'h"^-''^^'^'^"'^' dr^ 
Jo 

Jo 

= C12 

using the change of variable i] r]t~^^^. 

In order to extend Theorem El to second order differential operators other than 
Schrodinger operators, it needs to be generalized as follows. 

Theorem 19 Let k, : (0, 00) ^ C be a continuous function, let Cq, Cqo be positive 
constants and let aQ, Ooo be non-negative constants such that 

co'r'' < Re«:(0 < Cor° ^/ < ^ < 1, 
c-^r-" < Re«:(0 < if l<e<oo. 

Then the conclusion of Theorem^^ is still valid if we replace / IJ6|) by 
9h,u,d^) = exp{z^(x - u)/h - (x - uf/2hK{^)}. 

Proof We make obvious adaptations to the proof of Theorem El up to ()38|1 . which 
becomes 

/■CO 

F{h,s) = / /i~i/2^ReK(0)^/^exp{-C6(e//i-s)2/iRe/€(0} d^ 
Jo 

^0 

/oo 
h-y^c]i'C-/'exp{-ce{^/h-srhc-JC-} de 

Each of these integrals is estimated by the same method as in Theorem 1171 

9 Constructing the boundary pseudospectra 

When one examines the pseudo-eigenfunctions in several exactly soluble examples, 
ini UHl Unj; one sees that they do not conform to the above ideas. They are 
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strongly localized at one end of the interval in question, and decrease exponentially 
as one moves away from this end. 

In this section we develop the general theory of boundary pseudospectra for variable 
coefficient operators in the one- dimensional context. A partial extension to higher 
dimensions and manifolds is described in the next section. We assume that 

{Lhf){x) = -h^a{x)f"{x) - ihh{x)f'{x) + c{x)f{x) 

for X G [0,7]. The semiclassical principal symbol is 

a{u,i) = a{u)i^ + h{u)i + c{u). 

We will need the fact that the symbol can be analytically continued to complex ^, 
but only assume the coefficients of Lh, and therefore a, to be in u on [0,7]. 
Similar but weaker estimates can be proved if the coefficients are only C" for some 
n. We assume ellipticity, in other words that a(x) 7^ for all x G [0,7]. We start 
by ignoring the boundary conditions and looking for a pseudo-eigenfunction of the 
form 

f{s) = h-^/\{s)eM^{s)) (40) 

where 

m=—l 

We assume that x ^ C°°[0,7] satisfies xi^) = 1 if < s < 5/2 and = if 
s > 6] the constant 6 > must be small enough for the proof of Theorem l2T] to be 
valid. We put 

where 

w{^, v) = + ^ + + 

a{v) a{v) 4a(t')2 a{v) 

As before we take the branch of the square root which equals ^ + 6(0)/2a(0) at 
V = 0. However we now require Im(^) > 0, in order to ensure that f{s) decays 
rapidly as s increases. We have 

=i^s + ksy2 + 0{s^) 

for small s > as before. 

Lemma 20 Let F be a positive continuous function on [0, 6] and let G be a con- 
tinuous function on [0,6]. If m is a non-negative even integer then 



s'"G(s) exp{-h-\sF{s)} ds ~ ch"'+^ 

Jo 
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as h ^ 0+, where 

G{0)T{m + l) 



c 



m 



m+l 



In the following theorem we put {Qf){x) = xf{x) and (P/)(x) = —ihf'{x) as 
before. Although Q is self-adjoint on an obvious domain, we impose no boundary 
conditions on P, which is therefore not even symmetric. 

Theorem 21 // the coefficients of Lh are C°° and Im {^) > then for any positive 
integer n there exist functions f which depend on h, n, ^ such that 

Ihn 11/11 = c>0 (41) 

WQfW = 0{h) (42) 
\\Pf-m = 0{h) (43) 
\\L,f-a{0,Of\\ = 0{h-^') (44) 

as h ^ 0. 

Proof Let / be given by (^UJ. To prove 1)411) we write 

ll/ir = /^-^ /'x(5)^exp{2Re(^(.))}d. 
Jo 

= h'^ f G{s)exp{-h-hF{s)} ds 
Jo 

where 

F{s) = -2Re(V^_i(s))/s 
G{s) = x(s)'exp|2Re 

This is of the form treated by Lemma 1201 if 5 > is small enough to ensure that 
F{s) > for all s G [0,5]. 



To prove ()42|) we write 

||g/||2 = /i-i [\^G{s)exp{-h-hF{s)} ds 
Jo 

and apply Lemma QUI again. 

The proof of ()43|) uses Lemma EOl and the expansion 

Pf - ^f = 1^1+ 1^2 + /is 
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where 



/xi = -^h-'/'{i:'_,{s)-^0}x{s)exp{^p{s)} 
= -^h'/^(j2h^^|J'^{s)]x{s)exp{^|J{s)} 

\m=0 / 

/i3 = -ih^^'^x'is)exp{'ilj{s)}. 
The proof of ()44|1 follows in a similar way from the formula 

\m=ri+2 / 

We finally assume the boundary conditions 

uhf'{0) + wf{0) = (45) 

for some complex constants u, w, not both zero. We say that satisfies the 
exit condition at if Im (— 6(0)/a(0)) > 0. This language is motivated by the 
example discussed in in which Lh is the generator of a subMarkov diffusion on 
an interval. Given the exit condition at 0, we define the boundary semiclassical 
pseudospectrum at to be the set 

A = : < Im (0 < Im (-6(0)/a(0))}. (46) 

If ^1 G A and z = a{0,^i) then the other solution ^2 of <j{0,^) = z also lies in A. 
We have ^1 = ^2 if and only if 2; = c(0) — 6(0)^/4a(0). The set cr(0, A) is the region 
inside the parabola P = {a{0,t) : t G R}. 

Those familiar with UHl 1^ will observe the close relationship between the 
above and the winding number calculations there. At a qualitative level the given 
operator can be approximated near the end of the interval by the operator whose 
coefficients are frozen to the values which they have at the endpoint. Our theorem 
below provides quantitative flesh to this idea. It also provides the precise form of 
the relevant pseudo-eigenfunction, which is not easy to guess from the constant 
coefficient case. 

Theorem 22 Let Lh satisfy the exit condition at and let z lie inside the parabola 
P. Assuming z 7^ c(0) — 6(0)^/4a(0), let C,i, ^^2 ^ A denote the two distinct solutions 
ofcr{0,C,) = z. Given h > and n > 1, let fr be the boundary pseudo-eigenfunctions 
associated with h,n,^r as in ^Uj) and Theorem \21l and let 

f = {^<^ + w)h-{^ui,+w)h. (47) 
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Then f satisfies the boundary condition at and 

\\L,f~zf\\/\\f\\=Oih-^') (48) 

as h ^ 0. 

Proof The assumptions imply that fr satisfy the estimates of Theorem |^ from 
which fl^H|) follows. The proof that / satisfies (fi^ depends upon the identities 

/,(0) = /1-1/2 and /;(0) = 

10 Higher Dimensions 

The extension of the above ideas to higher dimensions needs more machinery. 
We are mainly interested in bounded regions in with smooth boundary, but 
since the proof of our main result depends upon choosing local coordinates around 
a boundary point rather carefully, we write down the argument in a manifold 
context. Let X be a smooth A^-dimensional manifold with boundary dX. Let X 
be provided with a volume measure dvol which has positive C°° density v{x) when 
restricted to any coordinate neighbourhood U. 

The natural differential d : ^"(X) -> C''~\T*X) is given within U by 

df{x) = idJix),...,dnfix)) 

and the adjoint operator d* : C"(TX) ^ C"-i(X) acts on a section g G C"(Tf/) 
by 

d*g{x) = —v{x) ^dj{v{x)g^ (x)). 

The differential operator Lh is determined by three coefficient functions, all as- 
sumed to be C°° and complex- valued on X; we write and T* in place of T^, (g) C 
and T* ® C below. We assume that a{x) : T* T^., b{x) G and c{x) G C for 
all X e X. Given h> and / G C°°{X) we then put 

{LJ){x) = d*{a{x) df{x)) - thb{x)-df{x) + c{x)f{x). 

Throughout this section a dot indicates the natural action of a covector on a tangent 
vector at some point of X. In the coordinate neighbourhood U the above formula 
may be written in the form 

{Lhf){x) = -h\-\x)d, {y{x)a^'\x)duf{x))-ihV{x)d,f{x) + c{x)f{x) 

using the usual summation convention, or in the form 

(L,/)(x) = -h\''\x)d,,kf{x) - ihV{h,x)d,f{x) + c(/i,x)/(x) (49) 
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where 



b>{h,x) = b'{x) + hb{{x) (50) 
c{h,x) = c{x) + hci{x) + h'^C2{x). (51) 

The set of all operators of the form ()49j) is invariant under changes of local coor- 
dinates. 

The symbol of is given by 

ahix, = h^a^'^{x)ijik + hV{h, x)Q + c{h, x) 

which is not an invariant expression: d* and both depend upon the choice of 
the density v. However the semi classical principal symbol 

= a^^\x)^,^k + Vix)^j+c{x) 
= a{x)^ ■ ^ + b{x) ■ ^ + c{x) 

is invariant under changes of local coordinates. 

The following theorem is a mult i- dimensional 'boundary' analogue of Theorem II II 
We expect that there is also a mult i- dimensional analogue of Theorem |221 We 
choose a point in dX, label it p, and choose a complex cotangent vector ^ at p. 
We require that Im (^) has zero dot product with any vector at p which is tangent 
to dX and positive dot product with any inward pointing vector at p. If ?7 is a 
coordinate neighbourhood around p we always assume that p is represented by the 
point e R^. 

Theorem 23 Let Lh be of the form where all of the coefficients in i5(J\} . 
i51]} are C°° functions on U . Let the complex cotangent vector ^ at & dX satisfy 
the conditions of the last paragraph. Then for every sufficiently small h > there 
exists fh G C°°{X) which vanishes outside a neighbourhood ofO whose radius is of 
order h^/'^, and satisfies 

hmJIMh = c>0 (52) 
\\LufH-<^,i)fhh = Oik'/') (53) 

as /i — > 0. 

Proof Let denote the set of x G R''^ for which x^ > and let R^ denote the 
set of X for which = 0. We choose local coordinates around such that 

f/ = {x e R+ : |x| < p} 



27 



and put 

dU={xER^ : \x\ < p} 

for some p > 0. We write x = {x',x^) where x' G R^^"^ and x^ G R. Our 
assumptions imply that C, = (^', where ^' is real and rj := Im (^^) > 0. 

Put a = 1/2 and 7 = (A^ + l)/4. Let (pi be a smooth function on R^~^ which 
equals 1 if \x'\ < 1 and if \x'\ > 2. Let 02 be a smooth function on [0, 00) which 
equals 1 if < x'^ < 1 and if x^ > 2. Let (f){x) = (f)i{x')(f)2{x'^)- Then the 
smooth function 

Mx) = /i-^e*""'«-X/i-"x) 
on U has support with the required property for all small enough h > 0. 
To prove ()52|) we observe that 

WfhWl ~ v{0)h-^^ [ 0i(/i-V)2d^-ix' re-2'^"''^"02(/i-"x^)'da;^ 
as /i ^ 0. 

The proof of depends upon writing 

7 

Lhfh - cr{0,Ofh = J2 9m 

m=l 

where 



91 = 




92 = 


-ih^-''-^a^'\x)^je'''~'^''(f)k{h-''x) 


93 = 


-ih^-"-^a^'''{x)^ke"'"~'^""(pj{'^~"^) 


9i = 


-h^-^''~^a''\x)e'^~'^-^^j,k{h-''x) 


95 = 


h-''{b^{h,x) - V{0)}^je'''~'^''(j){h-''x) 


96 = 


-ih^-'^-W{h,x)e'^~'^-''<Pj{h-''x) 


97 = 


h-''{c{h,x) - c(O)}e^'^~'«-^0(/i-°x). 



We estimate the norm of each of these as above, obtaining \\9r\\2 = 0{h°') for 
r = 1,5,7, ll^^lla = 0{h^-'') for r = 2, 3, 6 and ||^^||2 = 0{h^-^'') for r = 4. Given 
these estimates, the optimal value of a is 1/2. 
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We next impose boundary conditions of the form 

hu{x')n{x', 0) ■ df{x', 0) + w{x')f{x', 0) = 

for all x' G dU, where the complex-valued coefficients u, w are C°° on dU; we 
assume non-degeneracy of the boundary conditions at in the sense that m(O') 
and w{0') do not both vanish. The vector field n on [/ is supposed to be smooth 
and transversal in the sense that it has a non-zero inward pointing component at 
every point of dU. We use the associated fiow to construct local coordinates. In 
other words we choose local coordinates for which the boundary conditions can be 
written in the form 

h u{x')dNf{x', 0) + w{x')f{x', 0) = 0. (54) 

We say that the complex covector ^ at is admissible under the following condi- 
tions. We require that Im {C,) has positive dot product with any inward pointing 
vector at 0. We require that the same conditions hold for a second complex covec- 
tor ^ at 0. We require that z := o"(0, ^) = a{0, ^) and that ^ ■ t = ^ ■ t G R for any 
vector t which is tangent to dU at 0. In the local coordinates specified above we 
are fixing ^' = ^' G R^^^ and assuming that the two solutions C,^ and C,n of the 
quadratic equation 

a{0,{^',s)) = z 

in s G C both have positive imaginary parts. We say that satisfies the exit 
condition at if the set of admissible ^ is non-empty. 

Theorem 24 If G is an admissible covector and z = o"(0,^) then under the 
above conditions there exist Qh G C°°{U) satisfying the boundary conditions (|5^| ) 
and also 

{x G t/ : |x| < c'/i^/^} (55) 
c> (56) 

0{h^/^) (57) 
as h ^ 0. 

Proof We put 

gh{x) = a{x')fh{x) + a{x)fh{x) 

where 

/^(x) = /i-V'^"'^X/i-^/2a;) 
h{x) = /i-V""'«-^0(/i-i/2^). 



supp(5(/,) C 
limllc/ftlh = 

WLnQh-zghh = 
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In this equation 7 = (A^ + l)/4 and is as in the proof of Theorem 1221 Also 
^ = (^')^Af) and ^ = (^', ^at). The coefficients a, a are to be determined. Before 
continuing, we mention that in the case of Dirichlet boundary conditions we put 
dh = fh — fh, that is a(x') = —a{x') = 1; most of the calculations below are much 
simpler in this situation. 

It is immediate from the definition that 

hdNgh{x',0) = h-^a{x')^N + a{x')^N}e'''''^'"''Mh~^^^x'). 
It follows that Qh satisfies the boundary conditions provided 

iu{x'){a{x')C,N + a{x')^N} + w{x'){a{x') + a{x')} = 0. 
This is solved by putting 



a{x') = w{x')+iu{x)^ 



N 



a[x 



= —w{x') — iu{x')^N. 

Since C,N 7^ ^N, 01 and a cannot both vanish near 0'. 

The validity of (j55|) is immediate. To prove (j56p we note that 



00 



«(x')e'''"'«-X/i-^/2x) + «(x')e*^"'«-"0(/i-i/'x)| v{x) d^o; 
a(O')e'^"'«"0(/i-^/2x) + a(O')e''^"'«"0(/i-^/2x)|'t;(O) d^x 
<\>^{hr^l''x'f d^-^x' X 







d^ 



>2(/i-'/'x^)2 dx^ 



R^- 



s""? d.^ 



t;(0) 
> 



R^- 



d s 



a(0')e'^'^^" + a(0')e 



ds 



TV 



as — > 0. 

The proof of (jFfjl depends upon writing 

-^^ftfi-ft - zgh = A;i + ^2 + ^3 + K 
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where 



ki{x) 
k2{x) 

/C4(x) 



a{0'){LhUx)-zMx)} 
aiO'){Ljhix)~zhix)} 
Lf,[{aix')-aiO')}Mx)] 
L;,[{a(x')-«(0')}A(x)] 



and then estimating each term as before. 
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